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We study the magnetic interaction between two superconducting concentric mesoscopic Al loops, 
close to the superconducting/normal phase transition. The phase boundary is measured resistively 
for the two-loop structure as well as for a reference single loop. In both systems Little-Parks 
oscillations, periodic in field are observed in the critical temperature T c versus applied magnetic 
field H . In the Fourier spectrum of the T C {H) oscillations, a weak 'low frequency' response shows 
up, which can be attributed to the inner loop supercurrent magnetic coupling to the flux of the 
outer loop. The amplitude of this effect can be tuned by varying the applied transport current. 



I. INTRODUCTION 

In 1962, Little and ParksB measured a mesoscopic superconducting cylinder in an axial magnetic field. The su- 
perconducting critical temperature T c ($/$o) showed oscillations periodic in the normalized flux, with the period 
corresponding to the superconducting flux quantum <£>o = h/2e. These oscillations in T c (&/$o) are a straightforward 
consequence of the fluxoid quantization constraint, which was introduced by F. LondonEI. Fluxoid quantization can 
be easily understood by integrating the second Ginzburg-Landau (GL) equation for the supercurrenttau 
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along a closed contour. Here, j is the supercurrent density, v is the superfluid velocity, S is the phase of the complex 
order parameter \I> = 1^1 e tS , and A is the magnetic vector potential. Integration along an arbitrary closed contour 
yields the following equation: 

$' = $ + i <£m*v-dT=N$ , (2) 

where the fluxoid is quantized in units of $o = h/2e and <I> is the applied flux threading the area inside the 
contour. The integer number N is the phase winding number, or also called the fluxoid quantum number, counting 
the number of flux quanta $0 penetrating the enclosed area. When the applied flux <!> is not equal to integer times 
the flux quantum <E>o, a supercurrent j has to be generated in order to fulfill Eq. (^|). 

For a superconducting ring of radius r, made of wires of vanishing width (w = 0), a one-dimensional (ID) GL model 
can be used to describe the onset of superconductivity. The relative T c variations can be written as: 



T c0 r 2 V $ o/ 

where £(0) is the coherence length at zero temperature, T c q is the critical temperature in zero field , and the integer 
number N is chosen to maximize the critical temperature T c (&). For each fluxoid quantum number N, the critical 
temperature T c ($/$o) has a parabolic shape. The Little-Parks (LP) oscillations appear due to the transitions from 
the integer value N to N + 1, at a half integer value of &/&q. For <!>/<I>o = N, no supercurrent flows in the ring, 
while for &/&q = N + 1/2, the supercurrent reaches a maximal value and changes sign. The maximum normalized 
variation of the critical temperature AT c ($/$o)/T"cO is £ 2 (0)/4r 2 , as can be evaluated from Eq. To observe these 
quantization effects experimentally, structures with the radius of the loop of the order of the coherence length £(0) 
should be used. 

Since the pioneering work of Little: and Parks, single mesoscopic superconducting loops and cylinders has been 
largely studied. Recently, Zhu et aln studied the flux state in twcumagnetically coupled mesoscopic normal loops. 
The magnetic coupling of an array of normalQ and superconductingl] loops has also been studied. In those cases the 
loops are electrically isolated from each other and can only interact magnetically. 
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Our work will focus on the magnetic coupling of two concentric superconducting loops. The modification of the phase 
boundary T C {H) of the outer loop due to magnetic coupling with the inner loop will be studied. Magnetic coupling 
effects between two loops are potentially very important since for the two loops made from different materials, a new 
unusual effect of enhancing T c i(H) due to a higher T C 2(H) can be expected. 

The GL free energy for such system can be written as follows: 
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with T n the total free energy in the normal state, a and /3 the expansion coefficients, m* the mass of a Cooper pair 
and L and M the self- and mutual inductance. V is the volume of the loop, 4" is the order parameter, v is the velocity 
of the superfluid and I is the supercurrent. The indexes i and o refer to the inner and outer loop, respectively. The 
superfluid velocities of the two loops are determined from the fluxoid quantization constraint: 
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To solve this equation, the free energy must be minimized with respect to variations in ^>i, ty , Ii and I a . The free 
energy J- s contains in 
parameters ^ and ^ c 



energy T s contains in this case the term MIiI Q tx M \ ^>i\ 2 \ ^o\ 2 responsible for the mixing of the two individual order 



II. EXPERIMENTAL RESULTS AND DISCUSSION 




Sample A, B Sample C Sample R 

FIG. 1. (a) A schematic drawing of samples A, B and C, and sample R. Ii and I2 are the contacts for the transport current; 
Vi and V2 are the voltage probes, a; and a are the width of the inner, resp. outer loop, measured from middle to middle, 
and w is the wire width of the loop and of the connecting wires, (b) AFM image of the double loop (sample A, B and C) and 
of the single loop (sample R). 

We present the results of transport measurements carried out on two different types of mesoscopic Al structures 
(Fig. |l|). The first type of sample is composed of two concentric loops, with the outer loop being electrically connected 
to the experimental set-up in order to perform four-points resistance measurements. Three samples of this type have 
been studied: sample A and sample B with the same thickness and the same dimensions of the loops and sample C 
with a smaller thickness and slightly different dimensions of the loops. The reference sample (sample R) is analogous to 
the first structure, but without the inner loop. All the samples discussed, except sample C, are evaporated in the same 
run. All samples have been prepared by thermal evaporation of 99.9999% pure Al on a Si02 substrate. The patterns 
are defined using electron beam lithography on a bilayer of PMMA/PMMA-coPMMA resist before the deposition of 
an aluminum film with a thickness r=50 nm and r=28 nm for sample A, B and R and for sample C, respectively. 
After the evaporation, the lift-off was performed using dichloromethane. In Fig. [f]a, the geometry and the dimensions 
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of the different structures are shown. The thickness and the lateral dimensions of the samples have been characterized 
by the X-ray diffraction on a co-evaporated plane film and AFM (see Fig. [j]b), respectively. The wire width of the 
loop and of the connecting wires has been determined from SEM investigations. The superconducting/normal phase 
boundaries are obtained from transport measurements, carried out with a transport current it flowing through the 
outer loop. The phase boundary is measured holding the resistance at a fixed resistive criterion (we used the criterion 
Rn/2, with R n the resistance in the normal state). This is achieved using an electronic feedback circuit. Once a 
required temperature stability is obtained, the magnetic field is swept at a very slow rate (with a typical frequency of 
20 /iHz). To improve the signal-to-noise ratio a PAR 124A lock-in amplifier has been used, operating at 27.7 Hz. 

The width oLjbhe wires of the studied superconducting loops determines the parabolic background of the T C (H) 
phase boundary!!!: 
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where w is the width of the wires. Eq. (gj) also describes the T C (H) line for a superconducting thin film of thickness 
w subjected to a magnetic field parallel to the film plane. In a practical situation, the phase boundary of a supercon- 
ducting loop of finite wire width will show a parabolic background (see Eq. (S|)). The suppression of T c can be written 
as the sum of two components: an oscillatory term as described by Eq. (|3), and a monotonic term (Eq. (^J)). The 
coherence length £(0) of the samples can be determined from the parabolic background of the phase boundary, since 
the width of the wires if is a known parameter. A second method to evaluate the coherence length is to determine 
£(0) from the slope of T C (H) of a co-evaporated macroscopic reference film. A third method is based on the dirty limit 
of the GL theoryEl, where from the known value of the elastic mean free path I, the coherence length £(0) = 0.86yft^l 
is obtained, using the clean limit value £o=l-6 fJ-m for Al. The results of the three methods are summarized in table 
H The difference between the £(0) values calculated with the three different methods can be partially explained by 
the rather broad error margins on the wire width w. 
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FIG. 2. Calculated phase boundary T C (H) of the inner (solid line) and outer loop (dashed, dotted, and dashed-dotted line) 
for two periods of the inner loop, without magnetic interaction, using the dimensions of sample A and B. The phase boundary 
of the outer loop is shown for three different transport currents I t : lt=0 /-*A: dashed line, It— 0.3 /jA dotted line, and It —0.9 
(j,A dashed-dottedi4ip«. The coherence length £(0) = 103 nm is taken. The shift of the T C (H) curves with increasing It are 
estimated from reMB: I c = J c0 (T c0 - T c ) 3/2 ,with I c0 ps 550 /xA. 

In Fig. ||, the theoretical phase boundaries are shown for the two single loops when they are not coupled magnet- 
ically: one using the dimensions of the inner loop to calculate the phase boundary (solid line), and the other using 
the size of the outer loop (lowest dashed line), thus corresponding to the reference sample; this has been calculated 
from Eq. (|) using the dimensions of the loops summarized in table |. It should be mentioned that an increase 
along the vertical axis corresponds to a decreasing temperature. The amplitude of the LP oscillations is the larger 
the smaller the loop size (see Eq. (p5|)). The period of the T C {H) oscillations for the inner and outer loop is given by 
/^oAifi = 3>o/Si, an d HoAH = < f>o/S'o, respectively with Si — af and S a = al the enclosed areas (see Fig. ||). Points 
X and Y in Fig. are situated at two fixed applied magnetic field values on the T — H phase boundary for the outer 
loop. In the case of point X, the inner loop is in the normal state, thus the inner loop carries no supercurrcnt. In this 
situation, no flux is coupled to the outer loop. For point Y, on the contrary, the inner loop is in the superconducting 
state, and a supercurrent will flow in both loops in order to satisfy the respective fluxoid quantization constrains (Eq. 
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B 
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r (nm) 


50 


50 


50 


28 


w (nm) 


180 


180 


180 


140 


a (/Ltm) 


2 


2 
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1.9 


a, (/Ltm) 


- 


0.9 


0.9 


1 


T c0 (K) 


1.324 


1.326 


1.327 


1.363 


Ra (f2)=p/r (4.2 K) 


0.65 


0.72 


0.76 


1.36 


Zei (nm) 


12.3 


11 


10.6 


10.5 


£(0) (Eq. @) (nm) 


103 


105 


102 


117 


£(0) Ref. Sample (nm) 


128 


128 


128 




£(0) dirty limit (nm) 


114 


112 


120 


112 


U (pH) 




1.6 


1.6 


1.9 


L ( P H) 


4.5 


4.5 


4.5 


4.6 


Af (pH) 




22 


22 


21 



TABLE I. Material parameters for the measured samples 



). Under these conditions, due to the mutual inductance between the two current loops (see Eqs. (|j) and @), an 
influence of the fluxoid quantization in the inner loop on the measured T C (H) phase boundary of the outer loop is 
expected. 

To extend the flux interval for which the inner loop remains superconducting, a higher transport current It can be 
applied to the outer loop. The phase boundary T c ($) of the outer loop is schematically presented in Fig. || by the 
dotted (It=0.3 /iA) and the dashed-dotted (It=0.9 fiA) lines. We now follow through the fixed magnetic fields lines 
following the points X— >X" and Y— >Y" by increasing the transport current /(. The point on the phase boundary of 
the outer loop, with the same magnetic field value as X for lt=0, will cross the phase boundary of the inner loop 
in point X' by increasing I t . For a highest transport current, the inner loop will be superconducting for each point 
on the phase boundary of the outer loop (dashed-dotted line). The inner loop will be deeper in the superconducting 
state, following the shift from Y to Y" while increasing the transport current. As a result, an increase of It will not 
only broaden the interval in which the inner loop is superconducting, but also increase the supercurrent in the inner 
loop. 

Zhang et alB have calculated the self-flux for a typical mesoscopic ring. From these calculations, we can find the 
self- flux and the additional flux in the outer loop due to the presence of the inner loop. The self-flux for the inner and 
outer loop will be not larger than 0.4% of the applied flux for T/T c > 0.99. The additional flux Mli in the outer loop 
will be less than 4% in this temperature interval but can be higher than 20% for T/T c < 0.95. All the measurement 
presented in this paper are in the region T/T c > 0.99. A possibility to further increase the supercurrent in the inner 
loop Ii and thus the magnetic coupling between the two loops would be the use of a material with a higher T c for the 
inner loop. 

In Fig. ||, the phase boundaries of the samples A, B, C and R (□: sample A, : sample B, A: sample C, V : sample 
R) are shown for a ac transport current It =0.3 fiA rms after subtraction of the monotonic parabolic background due 
to the finite width of the strips (Eq. (6)). Further on, the flux $ = S /j,qII will refer to the flux threading the surface 
of the outer loop. The T c ($) part below 5 $o is not shown because the experimental data were rather noisy in the 
low field region for some measurements. The curves corresponding to different samples are arbitrarily s hifte d in Fig. 
|3|. The variation of the critical temperature with increasing transport current can be estimated from reiL3liJ for zero 
field: I c — I c o(TcO ~ T c ) 3 / 2 , with 7 c o «550 //A and I c o «240 ^A for sample A, B and R and for sample C, respectively. 
The critical temperatures of the different samples after extrapolation of the experimental results to lt=0 are given 
in table |[ For 7 t =0.3 /xA, the shift of the zero field critical temperature AT c q (used in Fig. ||) is 7 mK (12 mK for 
sample C), for 0.5 /iA: 9 mK and for 0.7 /iA: 12 mK. Within each oscillation period, a parabolic function is fitted 
through the data points in between the transition to a different fluxoid quantum number N a — > N a + 1 (Eq. (J|)). The 
value for the self- and mutual inductance^ for the different samples are summarized in table |. 

In Fig. H[ a smaller amplitude of the oscillations is observed for sample C, compared to the other samples. This 
smaller amplitude is to be expected because of a smaller thickness of the film, what would lead to a reduced mean 
free path and a smaller coherence length value for sample C (see Eq. (||) ) . But the coherence lengths are comparable 
for all four different samples (see table |) . 

The phase boundaries of the double loop sample A and for the reference sample R are also measured for different 
transport currents I t . The results for three different ac transport currents I t (0 : 0.3 fiA, A: 0.5 fj,A, XJ: 0.7 fiA 
rms) are shown in Fig. Ek and Fig. 0b for sample A and sample R, respectively. From these experimental data, clear 
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FIG. 3. Measured phase boundary T c ($) after subtraction of the parabolic background caused by the finite wire width of 
the loops. The phase boundary is plotted in normalized units of the flux $ = S ^oH threading the outer loop. Within each 
oscillation period, a parabolic function is fitted through the data points in between the transition to a different fluxoid quantum 
number N . The phase boundary of sample A (□), sample B (O)i sample C (A) and sample R (v) is compared for the same 
transport current It= 0.3/xA. 



differences between sample A and sample R are seen. First of all, the amplitude of the T c ($) oscillations is stronger 
in sample R, secondly the phase boundary of the single loop (Fig. 0b) matches very well with the fitted parabolic 
curves (solid curves), while the oscillations of sample A (Fig. [|a) are not parabolic at all, since the cusps in T C (<I>) 
are always rounded. 

The rounding of the cusps in T c ($), observed for sample A, are not reproduced in samples B and C (see Fig. 0). 
These rounded cusps cannot be attributed with all certainty to the magnetic interactions with the inner loop, but 
may be also relateid,to the presence of scattering imperfections in the outer loop of sample A, as has been shown 
theoretically in reic3. A possible source of such imperfections coukLb£ the variation of the strip width along the 
loops written by e-beam lithography. Using the micronet approacHLJ'Ej, the GL equation can also be solved when, 
imperfections are present in a loop. A single period of the phase boundary calculated from this micronet approachE3 
is presented in Fig. for different magnitudes of imperfections. Notice that the stronger the imperfection the more 
the oscillation amplitude is damped. n 

It is clear from our experiments and fromta that the influence of the inner loop on the measured phase boundary is 
quite small. To detect traces of the periodicity coming from the T c ($) oscillations of the inner loop, a fast Fourier 
transform analysis of the phase boundaries is carried out using 2 10 equally spaced points in an interval between 5 
and 24 $/4> for sample A, B and R and between -36 and 49 <&/$() for sample C, but with approximately the same 
number of data points. In order to remove most of the contribution arising from the fluxoid quantization in the outer 
loop, the parabolas fitted in Fig. E3, 0a and |]b (solid line) are subtracted from the experimental data, prior to taking 
the Fourier transform. The resulting spectra are presented in Fig. 0a for sample A, 0b for sample B, 0c for sample C 
and 0d for the single loop (sample R), for an external transport current of I t —0.3 [iA. The Fourier spectra of sample 
A and sample R for I t —0.5 fiA are shown in Fig. 0a and [lib, respectively. 

In Fig. |a, 0c and 0a, a peak is clearly seen at ($/$o) = 1 corresponding to the frequency of the LP oscillations 
for the outer loop; the second to the sixth harmonics of this base frequency are also seen. In the case of the sample R 
(Fig. 0d and ^>) and of sample B (Fig. 0b), the first harmonic is strongly reduced. This reduction is quite reasonable 
because of the subtraction procedure applied before the Fourier transformation. The lower scale for the Fourier 
spectrum for sample C (Fig. 0c) is probably due to a smaller amplitude of the oscillations in the phase boundary. For 
the double loop (sample A, B and C), supplementary peaks are clearly distinguished in between the harmonics (see 
inset in Fig. 0a, 0b, 0c and 0a), which are considerably weaker in the reference sample R. These peaks can be due to 
the coupling with the supercurrent in the inner loop. Comparing the insets of Fig. |6p and Fig. |fjd, we note that the 
supplementary peaks are substantially sharper and higher for sample B than for sample R. 

The period of the LP oscillations for the outer loop is noAH a = 4>o/5' =0.525 mT and /ioA_ff =0.550 mT for sample 
A, B and R and for sample C, respectively, corresponding to the first harmonic peak at ($/$o) _1 = L The periodicity 
of the T c ($) oscillations for the inner loop (samples A and B) is /ioAiJ^ = $o/5i=2.6 mT, which is approximately 
5 times larger than /ioAiJ Q , since the surface of the inner loop is approximately 5 times smaller than the one of the 
outer loop. This periodicity is in good agreement with the measurement on sample B (Fig. 0b), where four peaks are 
indeed observed between the first and the second harmonic, thus indicating a five times smaller frequency compared to 
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FIG. 4. Measured phase boundary T c ($) after subtraction of the parabolic background caused by the finite wire width of 
the loops. The phase boundary is plotted in normalized units of the flux $ = S ^oH threading the outer loop. Within each 
oscillation period, a parabolic function is fitted through the data points in between the transition to a different fluxoid quantum 
number N . The results for three different ac transport currents It (O 1 0-3 ^A, A: v : 0.7 //A) are presented (a) for 

sample A and (b) for sample R. 




FIG. 5. The calculated phase boundary is shown for one period using the micronet approach for a loop containing one 
scattering imperfectiontj. The different curves correspond to different magnitudes of the scattering imperfection: from zero 
magnitude for the upper curve to a large magnitude for the lowest curve. 

the frequency of the T c ($) oscillations of the outer loop. The T c ($) measurements on sample A show two pronounced 
peaks between each harmonic (Fig. ^a and|7]a). This suggests a period for the inner loop ^AHi ss 3/x Ai/ o which is 
in disagreement with the dimensions of sample A. For sample C, which has a slightly different size, the periodicity of 
the inner loop oscillations in T C (Q) is fioAHi = <£>o/Si=2.1 mT« AfioAH Q . This periodicity is in agreement with the 
measurements where three peaks are clearly observed before the first harmonic (see Fig. oc). 

The reason of the disagreement between the periodicity of the oscillations and the peaks in the Fourier spectrum 
of sample A can be a different effective surface. To calculate the periodicity of the oscillations of the inner and outer 
loop, the average size of the loops (through the middle of the wires) has been used. It is possible that we have to 
take a slightly smaller or larger effective surface for sample A. The dimensions of samples A and B may also be not 
exactly the same. But the effective surface has to be taken larger than the largest dimensions of the inner loop to 
obtain a value of n^AHi sa 2,jj, AH o . We therefore think that the smaller peaks might be hidden, and not clearly seen 
however in the Fourier spectrum. In that case, a value for the periodicity of 6 times the periodicity of the inner loop 
is more realistic for sample A. Coming back to Fig. 0, we can see that the interval (between 5 and 24 $/$o)j where a 
Fourier transform was performed for sample A, B and R, corresponds only to 3 or 4 periods /loAiJ;. The resolution 
of the Fourier spectrum in the low frequency regime will therefore be low, thus presenting an additional difficulty in 
interpreting the intermediate peaks in the Fourier spectra. The phase boundary T C (<I>) of Sample C is measured over 
a broader interval (between -36 to 49 <&/ <&o)- This interval corresponds to more than 20 periods figAHi, what results 
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FIG. 6. Fourier transform of the phase boundary after subtraction of the fitted parabolas within each oscillation period, for 
a transport current I t = 0.3 /j,A for (a) sample A, (b) sample B, (c) sample C and (d) sample R. The inset shows a zoom of the 
plot for the low frequency region. 



in a better resolution of the Fourier spectrum in the low frequency regime. 

In a first approximation, with the currents It and I independent from each other, the mutual inductance M can 
be evaluated from the amplitude of the additional peaks in the Fourier spectrum. The additional energy due to the 
mutual inductance in Eq. |j (MIiI ) has to be of the same order of magnitude as the amplitude of the oscillations 
ksATcoupUng due to coupling in the phase boundary of the outer loop. The value for the mutual inductance can then 
be evaluated with the formula 



M 



coupling 



(Iilo) 



(7) 



with Ub the Boltzmann constant. This gives M ss 17 pH, 13 pH_and 14 pH for sample A, B and C, respectively. 
The average of the supercurrents Ii and Iq is calculated from rc£m, for 7 t =0.3 jiA and with the criterion for T co at 
90% of i?„, what correspond to a typical temperature interval of 3 mK between this criterion and the 50% criterion 
used for the measurements of the phase boundary. It is quite clear that the supercurrents, and thus also the mutual 
inductance are strongly dependent of this chosen criterion. The values calculated from Eq. ^ are comparable to the 
calculated mutual inductances from table |. Due to the strong temperature dependent currents, only the order of 
magnitude of the mutual inductance can be evaluated from our measurements. 

For high (7 t =0.9 /iA for sample A and 1.0 /iA for sample C) transport current, the peaks between the harmonics 
completely disappear (not shown). This vanishing may be due to increasing shift of T c (<&) with the applied transport 
current. The highest current corresponds to the case of the upper curve (dashed-dotted line) in Fig. ||. Hence, for 
all flux values, the inner loop is in the superconducting state at T C (<I>) of the outer loop, and the phase boundaries 
of the two loops are not intersecting each other. There will be no sharp interruption of the supercurrent of the inner 
loop going from the superconducting to the normal state. Therefore, the shape of the phase boundary of the outer 
loop can be less sensitive to the presence of the inner loop. On the other hand, once the inner loop is deep in the 
superconducting state, a higher supercurrent in the inner loop would be present, and this higher supercurrent would 
require a higher coupling. If the transport current It is high enough, the inner loop will be always superconducting 
at T c ($) of the outer loop. A discontinuity in the measured phase boundary of the outer loop is expected when the 
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FIG. 7. Fourier transform of the phase boundary after subtraction of the fitted parabolas within each oscillation period, for 
a transport current It= 0.5 /iA for (a) sample A and (b) sample R. The inset shows a zoom of the plot for the low frequency 
region. 



inner loop changes from a fluxoid quantum number N% to JV,- ± 1. However, we could not see a discontinuity in the 
measured phase boundaries, corresponding to a sign reversal of the supercurrent in the inner loop. 



III. CONCLUSION 



We have measured the normal/superconducting phase boundary of a superconducting system consisting of two 
concentric mesoscopic loops, to study magnetic interactions between the two loops. The modification of the T c {<&) 
oscillations of the outer loop is seen in the Fourier spectrum of the T C (Q>) line due to the coupling between the outer 
and the inner loops. To interpret these observations, we have used two different models. The first model assumes 
the presence of scattering imperfections in the outer loop. This model cannot explain the observed evolution of the 
Fourier spectrum with the current, although it might be applicable for a fixed weak current. The second model 
explains the extra peaks in the Fourier spectrum by the magnetic coupling of the two loops. The systematic shift of 
the T c (<i>) phase boundary of the outer loop with the applied current It induces a well defined evolution of the Fourier 
Spectrum which was indeed found in our experiments. This evolution of the extra peaks in the Fourier spectrum 
with the applied current gives an experimental evidence for the presence of the magnetic interaction between the two 
superconducting loops. 

Future magnetic measurements on huge arrays of magnetically coupled loops, deeper in the superconducting state, 
could be helpful to reveal an enhanced magnetic coupling of both loops at lower temperatures. A inner loop made 
from a different superconductor with a higher critical temperature would certainly increase the magnetic coupling 
between the two loops. An enhanced critical field is expected in this case for the lower T c loop, at the expense of 
sharing the fluxoid quantization "burden" with a loop where superconductivity is stronger. 
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